% 读取Excel文件
data = readmatrix('./input/红外汇总-特征变量提取绘图.xlsx');  % 将'your_file.xlsx'替换为你的实际文件名
% 从Excel文件中读取特征变量的索引
featureIndices = readmatrix('./input/红外汇总-特征变量索引.xlsx');  % 将'your_file.xlsx'替换为你的实际文件名

% 提取数据
x = data(1, :);  % 第一行作为x轴数据
y = data(2, :);  % 第二行作为y轴数据

% 提取标注的特征变量序列
x_index = featureIndices(1, :);
y_index = featureIndices(2, :);

% 创建图形窗口
figure;

% 绘制曲线
plot(x, y, 'b-', 'LineWidth', 1.5);  % 使用蓝色实线绘制，线宽为1.5

% 可选：绘制数据点
hold on;
% 随机选择部分数据点进行标记（这里选择约20%的数据点）
numPoints = length(x);
numToMark = round(numPoints * 0.007);  % 选择20%的数据点
randIndices = randperm(numPoints, numToMark);  % 随机选择数据点索引

% 绘制带透明度的方块标记
% 假设 x 和 y 是你的数据
% randIndices 是随机索引
h = scatter(x_index, y_index, 64, 's', 'filled', 'MarkerEdgeColor', 'none', 'MarkerFaceColor', 'r');
set(h, 'MarkerFaceAlpha', 0.5);  % 设置标记透明度为50%



% 添加图形元素
% title('数据点图');  % 设置标题
xlabel('Wavenumber /cm^{-1}');  % 设置x轴标签
ylabel('Absorbance');           % 设置y轴标签
grid on;  % 添加网格

% 可选：添加图例
legend('Speetral curye', 'Selected wavenumber by ( ARS )', 'Location', 'best');
